Collective mechanics of embryogenesis: Formation of ventral furrow in Drosophila 
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We propose a 2D mechanical model of the ventral furrow formation in Drosophila that is based 
on undifferentiated epithelial cells of identical properties whose energy resides in their membrane. 
Depending on the relative tensions of the apical, basal, and lateral sides, the minimal-energy states 
of the embryo cross-section includes circular, elliptical, biconcave, and buckled furrow shapes. We 
discuss the possible shape transformation consistent with reported experimental observations, ar- 
guing that generic collective mechanics may play an important role in the embryonic development 
in Drosophila. 

PACS numbers: 87.17.Pq, 87.19. lx, 87.19. rd 



The single most important step in the embryonic de- 
velopment of all multicellular animals is gastrulation [l| . 
It establishes the future body plan of the organism and 
transforms a topologically spherical embryo into a topo- 
logical torus. In the beginning of this crucial process, a 
part of the convex shell-like single layer of cells called 
the blastula folds in, and the invaginated part of the 
embryo develops into the digestive system The ge- 
ometry of invagination differs quite considerably among 
the animals as does the shape of the embryo itself [H-ll]- 
One of the best-studied types of gastrulation is the for- 
mation of the ventral furrow in the fruit fly (Drosophila 
melanogaster) 0-[gJ- At the beginning of this process, 
the Drosophila embryo is an elongated ellipsoidal onc- 
cell-thick epithelium of about 6000 cells coating the yolk. 
During the first stage of the furrow formation, a length- 
wise invagination occurs as the cells in a stripe on ventral 
side deform and buckle inwards [4] (Fig. QJi,b). This pri- 
mary invagination globally alters the shape of the embryo 
and is one of the vital factors for determining the future 
head-to-tail body axis 

The role of mechanics in the early embryonic develop- 
ment is universally appreciated Q , and a range of models 
have been proposed to explain the formation of ventral 
furrow. To the best of our knowledge, all of them rec- 
ognize that the invaginating cells — the mesoderm - 
are different from the rest of the embryo 0, . By 
forcing the mesoderm cells into a keystone shape with a 
small outer (apical) face and a large inner (basal) face, 
the ventral section of the embryo can be driven rather 
naturally into an invaginated shape. This can be done by 
superposing active and passive deformations; the former 
are generated by an inherent preferred keystone shape of 
cells and the latter result from the elasticity of the ep- 
ithelium 

S H-E3- Alternatively, invagination can be re- 
produced by hydrodynamic motion of embryo immersed 
in a surrounding fluid and consisting of cells character- 
ized by surface energy [ill ]. In this theory, formation of 




FIG. 1: Drosophila ventral furrow (micrograph reproduced 
with permission from Ref. 3) is the lengthwise invagination 
of the embryo on ventral side (panel a). The cross-sections of 
the central part (reproduced with permission from Ref. [2j]) 
show the inward buckling of the initially circular epithelial 
wall (panel b). In our 2D model (panel c), the cross-section 
consists of a ring of N quadrilateral cells of area A c enclosing 
the yolk of area A y . All apical, basal, and lateral edges are 
characterized by line tensions r a ,rb, and Fi, respectively. 



ventral furrow is triggered by an increase of the apical 
membrane tension of the ventral cells. These sophisti- 
cated models have been elaborated both in the simplified 
2D variants describing the predominantly cylindrical cen- 
tral section of the embryo [H, 0, as well as in the full 
3D geometry [fj], [Io[ , and they also included the effect of 
the vitelline membrane which covers the embryo [{j- 

In this paper, we challenge the importance of cell dif- 
ferentiation as the main postulate of the existing theo- 
ries of ventral furrow formation in Drosophila. Instead 
we explore the primary invagination as a collective phe- 
nomenon in a system of cells of identical mechanical prop- 
erties. In our 2D model of the embryo cross-section that 
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is based on surface energy alone, each type of cell side is 
characterized by a specific line tension: The apical, the 
basal, and the lateral line tensions are all different from 
one another but each of them is the same in all cells. We 
compute the phase diagram of the model epithelium sub- 
ject to constraints associated with the incompressibility 
of the yolk and each cell to find that it does include the 
invaginated shapes even though the model is devoid of 
any cell specificity. A more detailed analysis shows that a 
thick enough epithelium containing a large enough num- 
ber of cells is needed for the invagination to take place. 

We take our starting geometry at the time before ven- 
tral furrow formation begins and the central part of the 
embryo is cylindrically symmetric. Experiments show 
that the size of each cell along the long axis of the em- 
bryo is almost constant during the process [45 and the cell 
content is incompressible. If the cells are approximated 
by prisms oriented lengthwise relative to the long axis of 
the embryo, this suggests that the area of the cell sides 
normal to the embryo long axis is conserved. In this case, 
the cross-section of the central part of the embryo cap- 
tures the only variable terms in the Hamiltonian based 
on surface energy. 

In our 2D model, the N cells in the embryo cross- 
section are represented by quadrilaterals arranged in a 
ring that encompasses a given amount of yolk (Fig. Q}:) . 
Both the cells and the yolk are assumed incompress- 
ible [H, HH so that the area of each quadrilateral cell A c 
is fixed and identical for all cells, and the area of the 
polygonally shaped yolk A y is fixed too. Apart from in- 
compressibility, we disregard any bulk properties of cells, 
yolk, and the surrounding medium. We attribute the 
energy of the system entirely to the three types of cell 
sides 12, and the key assumption of our model is 



the total energy [Eq. (1)] can be written in dimensionless 
form 



that the line tensions of lateral, basal, and apical edges 
are not the same. This seems plausible, because the ten- 
sion on the lateral sides is determined by the cell cortex 
tension and on cell-cell adhesion, whereas the tensions on 
the apical and basal sides are controlled by the cell cortex 
tension and surface energy (which depends on the neigh- 
boring medium, and the yolk and the vitelline membrane 
are functionally dissimilar). 

The energy of the model epithelium reads 
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r a ,rf,, and Ti are the line tension of the apical, basal, 
and lateral edges of cells, respectively, and L l al L l b , and 
L\ are the lengths of these respective edges in cell i. The 
sum goes over all cells and the value of each of the three 
line tensions is the same in all cells. Using the reduced 
apical and basal line tensions 



w = al a + /3h + k, 
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where l a = Y.i L a/ R v^b = J2i L l/ R v and k = 
l/2'^ ji L % l /R y are the reduced sums of all apical, basal, 
and lateral edge lengths, respectively. The radius of the 
yolk before the formation of ventral furrow R y = -J Ay/n 
sets the length scale of the problem. 

To describe the formation of ventral furrow as closely 
as possible, we fix the number of cells N to 50 and set 
A c /A y = 1/60 w 0.0167 as suggested by the micrographs 
of the embryo cross-section [8j; thus the total area oc- 
cupied by the epithelium is 5/6 or about 83 % of the 
yolk area. The only remaining dimensionless parameters 
of the model are the reduced apical and basal line ten- 
sions a and /3, respectively. By numerically minimizing 
the energy [Eq. ([3])] subject to the two area constraints 
A c = const, and A y — const., we find the exact equilib- 
rium shapes of the blastula. Minimization is done using 
the Surface Evolver package [3] • 

Our central result is the diagram of stable shapes 
shown in Fig.[5J The (en, /3)-plane can be partitioned into 
three regions according to the symmetry of the shapes. 
At large a and (3 the stable shapes are circular whereas 
at a < 2 and /3 not too large the cross-section of the 
embryo is characterized by C2/1 symmetry. The smaller 
the reduced basal line tension, the more elongated the 
shapes in this region and the more pronounced the shal- 
low symmetric invaginations on the long sides. The third 
region of the diagram (shaded part of Fig. [2]) contains 
buckled shapes marked by a single invagination. Some of 
these are globally round just like the cross-section of the 
Drosophila embryo. For example, the (a = 2.2, (3 = 0.3) 
shape is quite similar to the experimental images of the 
embryo right after the primary invagination [i|,[l5j]. Many 
of our theoretical invaginated shapes are characterized by 
a very localized symmetry-breaking segment of the em- 
bryo wall. 

An intriguing feature of the shape diagram is the prox- 
imity of the circular and the invaginated shapes, which 
suggests that a direct transformation generating the in- 
vagination is possible. Depending on the starting state, 
one can think of many admissible pathways but a diago- 
nal cut involving a simultaneous decrease of both a and f3 
seems the simplest and thus the most appealing. Imagine 
starting with the (a = 2.4, f3 — 0.4) circular shape: De- 
creasing both reduced line tensions by about 13 % will 
transform it into a buckled shape as illustrated in the 
inset of Fig. [5J thereby reproducing the early stage of 
gastrulation in Drosophila. 

The reason why we choose to illustrate the ventral fur- 
row formation by a transformation involving a propor- 
tional decrease of the two reduced line tensions is rather 
transparent. The simplest way of materializing this par- 
ticular trajectory is to increase the lateral line tension 
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1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 
reduced apical tension a 

FIG. 2: Shape diagram of the embryo cross-section in the (a, /3)-plane for N = 50 and A c /A y = 1/60. The invaginated shapes 
qualitatively reminiscent of the Drosophila embryo after ventral furrow formation are located in the shaded region, and those at 
a around 2.2 and B no larger than 0.4 are quantitatively very similar to the experimental micrographs of embryo after primary 
invagination [l5j]. The rest of the diagram contains symmetrical elongated shapes (bottom left) and circular shapes (top 
right). A possible transformation of a circular shape into an invaginated one by proportionally decreasing the two reduced 
tensions a and 8 (A-E) is elaborated in the inset. 



while keeping the apical and the basal line tension un- 
changed. Thus the furrow formation can be induced by 
a variation of a single control parameter, and the varia- 
tion itself can be quite modest. A possible microscopic 
mechanism of increasing the lateral line tension may be 
a decrease of the cell-cell adhesion strength at fixed cor- 
tex tension. To obtain an even more buckled shape such 
as the (a — 2.0, f3 — 0.2) shape, a different mechanism 
where a and /3 do not change proportionally is needed. 

After having established that the proposed model pre- 
dicts shapes remarkably similar of those seen in the 
Drosophila embryo, we would like to understand its work- 
ings in more detail. Firstly, we vary the number of cells 
in the cross-section to see whether the discreteness of the 
model is essential. At fixed yolk and epithelium area, 
we halve N (thus setting A c /A y to 1/30 » 0.0333) and 
recompute the shape diagram. All other things being 
equal, this would halve the total energy associated with 
the lateral sides so it is not surprising that the shape di- 
agram for N — 25 can be mapped onto that for N = 50 
such that a given (N = 25,a25,/?25) shape corresponds 
to a (N — 50,cv5o = 2a 25 ,/3 5 o = 2/325) shape. The only 
exception are the furrow shapes absent in the N = 25 
diagram, which implies that a sufficiently large number 
of cells is needed to stabilize the buckled shapes. This 
conclusion is further supported by the N — 100 diagram, 



which does contain the furrow shapes and nicely agrees 
with the N = 50 results after the mapping analogous to 
that mentioned above. In turn, it suggests that there 
may exist a continuum limit of our model similar to that 
introduced in Ref. }13| . 

The circular, elliptical, and biconcave shapes shown in 
Fig. [2] closely resemble the shape of 2D lipid vesicles, in- 
finitely thin unstretchable but flexible loops [l6(. This 
is illustrated in more detail in Fig. [3] where four selected 
non-furrow shapes are replotted along with the contours 
of the matching 2D vesicles. The vesicle contours are 




a = 0.55 a = 0.70 a = 0.90 a = 0.95 

FIG. 3: Comparison of selected biconcave and elliptical model 
embryo cross-sections from Fig. [2] and the corresponding 2D 
vesicle contours, which are labeled by the reduced area a 
and closely match the neutral line of the epithelium. The 
epithelium cross-sections from left to right correspond to 
(a = 1.2,(5 = 0.2), (a = 1.4, 8 = 0.3), (a = 1.8, 8 = 0.5), 
and (a = 2.0, B = 0.5). In the part of Fig. [2] occupied by the 
vesicle-like shapes, the reduced area is constant along lines 
described approximately by a/2 + 8 = const. 



4 



characterized by the so-called reduced area a = 4irA/L 2 
(where A and L are the area and the perimeter of the 
vesicle, respectively) , and they were computed using Sur- 
face Evolver 14j. The very good agreement leads us to 
conjecture that the buckled shapes (which do not have a 
2D vesicle counterpart) must be stabilized by the finite 
thickness of the epithelium. 

To verify whether this is true, we varied the area of the 
cell while keeping the number of cells and the yolk area 
fixed. For N = 50, we computed the shape diagram for 
A c /Ay = 0.0333,0.025,0.0167,0.0133, and 0.00833. The 
results show that as the thickness of the epithelium is 
decreased, the region of the phase space where the buck- 
led shapes are stable gradually shrinks. For N = 50, the 
buckled shapes are present only in model embryo cross- 
sections with A c /Ay > 0.0133. In this sense, gastrulation 
in the elongated Drosophila embryo is very different from 
that in the spherical sea urchin embryo which can be 
qualitatively reproduced by a theory based on the elas- 
ticity of a thin fluid shell fl7j |. 

There are a few features of the ventral furrow forma- 
tion that are not captured by our model, the most im- 
portant one being the approach and the closing of the 
furrow (final stage in Fig. |T]b) . These processes can be 
phenomenologically included by attraction and adhesion 
of apical faces of cells on the opposite sides of the fur- 
row. In addition, the scenario proposed above involves 
equilibrium states of the epithelium and disregards the 
hydrodynamics of the furrow formation, which could be 
introduced like in Ref. 11]. Elaborating our model to 



account for these aspects of the process seems feasible, 
but simplifying it in any way without compromising its 
predictions is likely much more difficult. 

The main conclusion of this work is that from the me- 
chanical perspective, neither cell differentiation }18j nor 
the ensuing inhomogeneous active or passive deforma- 
tions are needed for the primary invagination in the ven- 
tral furrow formation. This is consistent with the at- 
tenuated furrow seen in mutants with suppressed meso- 
derm differentiation, which involves fewer cells and is not 
as deep as in normal embryos [191 ] ■ Moreover, our re- 
sults show that the buckled yet overall round shape is 
stable in absence of the vitelline membrane, suggesting 
that the ventral furrow in Drosophila can be shaped by 
the generic, spatially uniform, and collective mechanism 
proposed here. Its role relative to the localized morpho- 
genetic processes driven by cell differentiation is hard to 
assess, but the two are probably active simultaneously. 



If embryogenesis were engineered, differentiation would 
probably harness and steer the thrust provided by the 
buckling of the thick epithelium. 
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